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Abstract 

The heavy ion cross section for continuum e + e~ pair production has been calculated to all 
orders in Za. The formula resulting from an exact solution of the semiclassical Dirac equation in 
the ultrarelativistic limit is evaluated numerically. An energy dependent spatial cutoff of the heavy 
ion potential is utilized, leading to an exact formula agreeing with the known perturbative formula 
in the ultrarelativistic, perturbative limit. Cross sections and sample momentum distributions are 
evaluated for heavy ion beams at SPS, RHIC, and LHC energies. e + e~ pair production probabilities 
are found to be reduced from perturbation theory with increasing charge of the colliding heavy 
ions and for all energy and momentum regions investigated. 
PACS: 25.75.-q, 34.90.+q 
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I. INTRODUCTION 



The calculation of heavy ion induced continuum e + e~ pair production to all orders in 
Za is of continuing interest because up to now it has only been carried out for total cross 
sections and in a limiting approximation ]], Q, 0- Recent progress on this topic began with 
the realization that in an appropriate gauge |4|, the electromagnetic field of a relativistic 
heavy ion is to a very good approximation a delta function in the direction of motion of 
the heavv ion times the two dimensional solution of Maxwell's equations in the transverse 
direction |5]. This realization led to an exact solution of the appropriate Dirac equation 
for excitation of bound-electron positron pairs and a predicted reduction from perturbation 
theory of a little less than 10% for Au + Au at RHIC [6|. This reduction can be identified 
as a Coulomb correction to bound-electron pair production. 

It soon followed that an analytical solution of the Dirac equation was obtained inde- 
pendently and practically simultaneously by two different collaborations 0, ^, E| for the 
analagous case of continuum e + e~ pair production induced by the corresponding counter- 
moving delta function potentials produced by ultrarelativistic heavy ions in a collider such 
as RHIC. An extended discussion and reanalyis of this solution, with comments on early 
parallel work in the literature, shortly followed [lo|. Baltz and McLerranj^] noted the ap- 
parent agreement of the obtained amplitude with that of perturbation theory even for large 
Z. Segev and Wells|9] further noted the perturbative scaling with Z\Z\ seen in CERN SPS 
data 11]. These data were obtained from reactions of 160 GeV/nucleon Pb ions on C, Al, 
Pa, and Au targets as well as from 200 Gev/nucleon S ions on the same C, Al, Pa, and 
Au targets. The group presenting the CERN data, Vane et al., stated their findings in 
summary: "Cross sections scale as the product of the squares of the projectile and target 
nuclear charges." On the other hand, it is well known that photoproduction of e + e~ pairs 
on a heavy target shows a negative (Coulomb) correction proportional to Z 2 that is well 
described by the Bethe-Maximon theory [l^. 

Several authors subsequently argued [l|, Q, [J that a correct regularization of the exact 
Dirac equation amplitude should lead to a reduction of the total cross section for pair 
production from perturbation theory, the so called Coulomb corrections. The first anal- 
ysis was done in a Weizsacker- Williams approximation [|. Subsequently, Lee and Milstein 
computed 2], [J the total cross section for e + e~ pairs using approximations to the exact ampli- 



tude that led to a higher order correction to the well known Landau-Lifshitz expression [13J. 
In a previous paper [14] I have tried to explicate the Lee and Milstein approximate results 
and argued their qualitative correctness. 

In the present paper I undertake the full numerical calculation of electromagnetically 
induced ultrarelativistic heavy ion electron-positron pair production. I utilize a cross section 
formula derived from the exact solution of the Dirac equation with an appropriate energy 
dependent cutoff of the transversely eikonalized potential employed. 



II. CROSS SECTIONS WITH HIGHER ORDER COULOMB CORRECTIONS 

For production of continuum pairs in an ultrarelativistic heavy ion reaction one may 
work in a frame of two countermoving heavy ions with the same relativistic 7, and the 
electromagnetic interaction arising from them goes to the limit of two 5 function potentials 

V(p, z, t) = 8{z -t)(l- a z )hr{p) + 5{z + t)(l + a z )A+(p) (1) 

where 

iP/2? 



A±(p) = -Zamie"££. (2) 



The previously derived semiclassical amplitude for electron-positron pair production 0, 
B Q lo| written in the notation of Lee and MilsteinQ] takes the form 

M(p,q) = J -|^exp[zk-b]>l(k)F B (k)F A (q ± + p ± -k) (3) 

where p and q are the four-momenta of the produced electron and positron respectively, 
P± = Po ± Pz, q± = qo ± Pzi 7± = 7o ± Izi ol = 7o7; k is an intermediate transverse 
momentum transfer from the ion to be integrated over, 



w,, . — / \ a (k~P±) +lom 

M{k) = u{p) — —7-u{-q) 

~P+Q- — (k — pj_) z — m z + le 



taw "'r^ +T, ;. ^(-A (4) 

—P-q+ — (k — q±) — m 2 + ie 

and the effect of the potential Eq. (1-2) is contained in integrals, Fg and Fa, over the 
transverse spatial coordinates 0, US, Q|, 

F(k) = J d 2 p exp[-ik- p]{exp[-i2Za\np] - 1}. (5) 



F(k) has to be regularized or cut off at large p. How it is regularized is the key to under- 
standing Coulomb corrections. 

Mthough, a 8 has bee, pointed outQ, the derived exact se—ca. Dhac a.plitude 
is not simply the exact amplitude for the excitation of a particular (correlated) electron- 
positron pair, there are observables, such as the total pair production cross section, that 
can be constructed straightforwardly from this derived amplitude 
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18, 



3- 



This point 



has a long history of discussion in the literature |20L l2lL l22l |23(. The exact amplitude for a 
correlated electron-positron pair will not be treated here. The point is that exact solution 
of the semi-classical Dirac equation may be used to compute the inclusive average number 
of pairs — not an exclusive amplitude for a particular pair. Calculating the exact exclusive 
amplitude to all orders in Za is not easily tractable due to need for Feynman propagators [15]. 
The possibility of solutions of the semi-classical Dirac equation is connected to the retarded 
propagators involved. In this paper we do not consider the exclusive (Feynman propagator) 
amplitude at all. We concentrate on observables that can be constructed from the above 
amplitude and investigate the Coulomb corrections contained in them. 

In a previous article IjJ I have discussed these matters in more detail. There the uncor- 
related cross section expressions for da(p), da(q), and ot were presented, 

do{p) = I I ^|A1(k)| 2 |F.4(qx + Px - k)|VB(k)| 2 , (6) 

m 2 d?pd?q f d 2 k 



°r = J W^/ ^I*WIV.«U + Px - k)| (8) 
da(p) is the cross section for an electron of momentum (p) where the state of the positron 
is unspecified. Likewise da(q) is the cross section for a positron of momentum (q) with the 
state of the electron unspecified. Note that o T corresponds to a peculiar type of inclusive 
cross section which we should call the "number weighted total cross section" , 



[ d 2 bN= j d 2 bjr n p n (b), (9) 

3 J n=l 



in contrast to the usual definition of an inclusive total cross section 07 for pair production, 

o!= j d 2 bjrp n (b). (10) 



n=l 



In Eqs. (6-8), it is assumed that the sums have been taken over the electron and positron 



polarizations in |A^(k)| 2 . Taking traces with the aid of the computer program FORM 24] 
one obtains 

\M(k)\ 2 = 2p + g_[(k-p ± ) 2 + m 2 ] + 2p^q + [(k - q L ) 2 + m 2 ] 



m 2 [p + q^ + (k — p_i_) 2 + m 2 } 2 m 2 [p_g + + (k — q_i_) 2 + m 2 ] 2 

4[k • p±q+q- + k ■ q±p+p_ - 2k • p ± k ■ q ± + k 2 (p ± • q ± - m 2 ) - p+P-q+q-] 



m 2 [p + q_ + (k — p±) 2 + m 2 }[p_q + + (k — q_±) 2 + m 2 ] 
This expression exhibits the expected property that |.M(k)| 2 vanishes as k goes to zero; the 
positive squares of the direct and crossed amplitudes (terms one and two) are cancelled by 
the negative product of direct and crossed amplitudes of term three. These background terms 
can be subtracted off analytically from the expression for |.M(k)| 2 to obtain an expression 
exhibiting only terms dependent on k in the numerators: 

2D 2 r ]ll -2A ll {2D l f3 1 + (3 2 ) , 2D 2 V22 - 2A 22 {2D 2 {3 2 + (3 2 ) 



\M(k)\ 2 = 7, + 



m 2 D 2 (D 1 + A) 2 m 2 D 2 (D 2 + (3 2 ) 2 

DiD 2 r]i 2 - A 12 {D 2 p 1 + Dxjh + (u) 
m 2 D l D 2 (D 1 + p 1 )(D 2 + (3 2 ) [ } 



where 



An = p+q-(p± 2 + m 2 ) A 22 = P-q+(q± 2 + m 2 ) 

A 12 = -{p ± 2 + m 2 )(q ± 2 + m 2 ) 

D\ = p+q- + Pa 2 + m 2 D 2 = p_q + + q ± 2 + m 2 

fa. = -2k • p ± + k 2 f3 2 = -2k • q ± + k 2 

rj n = p+q~(3i r/22 = P-q+fc 

T]i2 = k ■ p±q+q- + k • q±p + p- - 2k • p ± k ■ q ± + A; 2 (p ± ■ q ± - m 2 ). (13) 

Every term in the numerators now has at least a linear dependence on k. This subtraction 
turned out to be necessary to limit roundoff error in calculations at the highest beam energies 
such as LHC. 

^, [n]] apart from a 



If one merely regularizes the integral itself at large p one obtains 
trivial phase 

. AnaZ 

F (V = ( 14 ) 
Then all the higher order Za effects in M(p, q) are contained only in the phase of the 
denominator of Eq. (14). Then, since the cross sections Eq. (6-8) go as |F(k)| 2 the phase 



falls out of the problem and it directly follows that calculable observables are equal to 
perturbative results. However, in this approach a lower k cutoff at some u/j has to be 
put in by hand to obtain dependence on the beam energy and to agree with the known 
perturbative result in that limit. 

Our present strategy is to apply a spatial cutoff to the transverse potential x{p) (which 
has been up to now set to 2Zahip) in order to obtain an expression consistent with the 
perturbation theory formula 13 m the ultrarelativistic limit. Instead of regularizing the 
transverse integral itself Eq. (5) and letting the cutoff radius go to infinity as was originally 
done Baa |10|, we will rather apply an appropriate physical cutoff to the interaction 
potential. In the Weizsacker- Williams or equivalent photon treatment of electromagnetic 
interactions the effect of the potential is cut off at impact parameter b ~ where 7 is 
the relativistic boost of the ion producing the photon and u is the energy of the photon. If 



and V(r) is cut off in such a physically motivated way, then 



X (p) = / dzV(y/? + p 2 ) (15) 

J — 00 

G 



where 



V(r) = - Za ^M-^ A ,B/l] (16) 



p+ + q+ P- + q- f-.^ 
= — ; uj b = (17) 



with uja the energy of the photon from ion A moving in the positive z direction and u>b the 
energy of the photon from ion B moving in the negative z direction. Note that we work 
in a different gauge than that used to obtain the original perturbation theory formula, and 
thus our potential picture is somewhat different. The transverse potential will be smoothly 
cut off at a distance where the the longitudinal potential delta function approximation is no 
longer valid. 

The integral Eq. (15) can be carried out to obtain 

X {p) = -2ZotK Q {pu AtB h), (18) 

and Eq. (5) is modified to 

F AB (k) = 2tt| dppJ (kp){exp[2tZ AB aK (puj AtB / 1 )} - 1}. (19) 



Fji(k) and Fs(k) are functions of virtual photon uoa and ujb respectively. The modified 
Bessel function Ko(pu/j) = — ln(p) plus constants for small p and cuts off exponentially at 
p ~ j/u- This is the physical cutoff to the transverse potential. 

One may define £ = kp and rewrite Eq.(19) in terms of a normalized integral I A ,B{lk/uS) 

F a ,b{\l) = A7lZ £ Ba i A Alk/cu) (20) 

where 

I A .B(lk/uj) = — ^ / dZZJ (Z){exp[2iZ AtB aK (fr/>yk)] - 1}. (21) 

It is now clear that Fa,b is a function of ArcZA^B/k 2 times a function of (jk/u). The limit 
as Z — ► of I^g^k/u) is analytically soluble 

'U^M = 1+ J 7 , V . (22) 
and one has ^^(k), the familiar perturbation theory form 

F ^ (k) " FT^W (23) 



As I have shown in a previous paper one might use some other physical cutoff and still 
obtain the Lee-Milstein Coulomb correction as long as one was expanding the Coulomb cross 
section correction only to lowest order in k 2 . However such an alternate physical cutoff would 
not lead to this correct perturbation theory form for F^ B (k) and would lead to modified 
results for the Coulomb corrections in a full integration over k. 

Fig. (1) displays the results of numerical calculation of ^(k^/uj)^ for Z = 82 and in the 
perturbative limit. Note that the upper cutoff of p at 7/0; has the effect of regularizing F(k) 
at small k. -F(k) goes to the constant Att^/uj 2 as k goes to zero in the perturbative case; 
it goes to a reduced constant value as k goes to zero for Z = 82. The form of the original 
solution Eq. (14) 

*W = ( 24 ) 
is simply wrong because it is unphysical. Since it lacks a proper physical cutoff in p, it not 
only blows up at k = 0, but it also fails to exhibit the correct reduction in magnitude that 
occurs when k'y/u is not too large. 

It is clear from Figure 1 that for large Z = 82 Coulomb corrections reduce |F(k)| 2 from 
the perturbative result for kj/u << 100. Only for k > ~ 100 does the magnitude of 
F(k) go over into the perturbative result. 




FIG. 1: The sold curve is the normalized integral squared for Z=82. The dashed line is the 
corresponding perturbation theory result. 

III. CALCULATIONS:NUMERICAL TECHNIQUE AND RESULTS 

The expression for the total cross section, Eq. (8), involves an eight dimensional integral 
over the positron and electron momenta as well as the virtual photon transverse momen- 
tum. This integral reduces to seven dimensions in the usual way by symmetry, e.g. let the 
positron transverse momentum define the x-axis. The usual method of evaluation, e.g. in 
perturbation theory, is via Monte Carlo. However, I have chosen to do the seven dimen- 
sional integral directly on meshes uniform on a logarithmic scale in each radial momentum 
dimension. It was possible to carry the calculation out without using Monte Carlo because 
the integrand is very smooth and smoothly goes to zero at both high end and low end of 
the momentum ranges. No artificial cutoffs were applied. 
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Calculations labeled exact and perturbative differ only in the expressions used for Fa,bQz). 
The analytical expression Eq. (23) is used for perturbative calculations. The exact calcu- 
lations makes use the expression Eq. (19), which must be evaluated numerically, but only 
once for each Z AB °f interest. 

Results of the numerical calculations will be compared with previously derived closed 
formulas for total cross sections. It is useful here to review those formulas. The Racah 
formula for the total e+e~ cross section in perturbation theory isH 



[Z 1 a) 2 (Z 2 a) 



7T77r 



28 
27' 



C 3 



178 £2 370 + 7tt 2 £ 



116 13tt 2 7 



where 



C = log 



27 



Pl-P2 



27 



9 



54 



9 



C(3) 



[lo g 2(2 7 2 -l)], 



(25) 



(26) 



M t M 2 . 

the relativistic 7 is that of each colliding ion in an equal and opposite ion velocity frame, 
and the Riemann zeta function 



00 1 

c(3) = E 



1.2020569. 



n=l 



n" 



(27) 



13], 



The log (7 2 ) term is the same as the original Landau-Lifshitz formula[13j, but the other ad- 
ditional terms are an improvement that allows this very early formula to attain a remarkable 
degree of accuracy as demonstrated by comparison with recent Monte Carlo evaluations. 



The Lee-Milstein formula^ includes higher order aZ effects in addition to the log 3 ( 7 2 ) 



term, 



{Z l a) 2 {Z 2 a) 2 2% 



27 



log 3 ( 7 2 ) - 3[f(Z A a) + f{Z B a)} log 2 ( 7 2 



where 



+%f{Z A a)f{Z B a)\og{Y 



f{Z AtB a) = Z 2 A B a 2 E 



(28) 



(29) 



n= -n(n 2 + Z\ B a 2 \ 

The dominant (negative) Coulomb correction in this formula is the log 2 (7 2 ) term, which was 
originally obtained by Ivanov, Schiller, and Serbo[l| with the Weizsacker- Williams approx- 
imation. The last (positive) log(7 2 ) term in the formula can be thought of as representing 
the Coulomb correction corresponding to multiple photon emission of both ions |28( and as 
we show below is relatively small. 



TABLE I: Computer calculations compared with analytical formula results. 7 is defined for one of 
the ions in the frame of equal magnitude and opposite direction velocities. Total cross sections are 
expressed in barns. The positive contribution of multiple photon emission from both ions to the 
overall difference betwen exact and perturbative results is shown in parentheses. 





Exact 


Perturbative 


Difference 


Pb + Au 


Computer Evaluation 


2670 


3720 


-1050 (+80) 


7 = 9.2 


Racah Formula 




3470 






Lee-Milstein 


3050 


5120 


-2070 (+160) 


S + Au 


Computer Evaluation 


119.7 


141.6 


-21.9 (+0.15) 


7 = 9.2 


Racah Formula 




132.0 






Lee-Milstein 


152.0 


195.0 


-43.0 (+0.30) 


Pb + Pb 


Computer Evaluation 


3210 


4500 


-1290 (+100) 


7 = 10 


Racah Formula 




4210 






Hencken, Trautmann, Baur 




4210 






Lee-Milstein 


3690 


6160 


-2470 (+190) 


Au + Au 


Computer Evaluation 


28,600 


34,600 


-6,000 (+220) 


7 = 100 


Racah Formula 




34,200 






Hencken, Trautmann, Baur 




34,000 






Lee-Milstein 


34,100 


42,500 


-8,400 (+290) 


Pb + Pb 


Computer Evaluation 


201,000 


227,000 


-26,000 (+600) 


7 = 2960 


Racah Formula 




226,000 






Lee-Milstein 


226,000 


258,000 


-32,000 (+700) 



Table I shows the results of numerical calculations. The present perturbative computer 
calculations are in good agreement with the Racah formula at RHIC and LHC energies as 
expected, and with the published Monte Carlo RHIC calculations of Hencken, Trautmann, 



and Baur 



29 



[3^ . At SPS energies the present perturbative computer calculation results 
are a bit higher (7%) than the Racah formula and the Hencken, Trautmann, and Baur 
calculation, perhaps indicating divergence in those results from the ultrareletivistic limit 
of the present treatment. The full numerical evaluation of the exact semi-classical total 
cross section for e + e~ production with gold or lead ions shows reductions from perturbation 
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theory of 28% (SPS), 17% (RHIC), and 11%(LHC). Clearly with increasing beam energy 
(and a larger value for the spatial cutoff of the transverse integral in the formula) higher 
order corrections to perturbation theory are relatively smaller. The S + Au calculation 
at SPS energy shows an expected smaller reduction from perturbation theory (15%) than 
the 28% reduction of Pb + Au at the same energy. The Lee-Milstein higher order overall 
correction to perturbation theory (difference column) is negative but somewhat larger than 
the difference evaluated here numerically. The small positive contribution of multiple photon 
emission from both ions to the overall negative Coulomb correction is shown in parentheses 
in the difference column. Because of the way the numerical calculations were organized it 
was straightforward to extract this contribution from the exact computer evaluation. Again 
the Lee-Milstein formula overestimates this small positive contribution, especially for the 
SPS case. 

There is the question of whether Coulomb corrections might become vanishingly small in 
some momentum regions. Let us take the Au + Au RHIC case as an example. If one looks 
at the uncorrelated positron cross section cross section (Eq. (7)) as a function of momentum 
then one finds that throughout the transverse and longitudinal momentum space of the final 
positron, the smallest reduction from perturbation theory is 12.5% and the largest reduction 
is 25% in comparison to the mean or integrated total cross section reduction of 17% of the 
table. Thus the argument given in Ref. |l| that Coulomb corrections contribute mostly for 
qx = m e , but should disappear for larger and smaller q^ is not verified. 

Figure 2 shows the transverse momentum spectrum integrated over all longitudidal mo- 
menta. The overall contribution does peak at about q^ = m e . However, Coulomb correc- 
tions persist to the highest and lowest values of q_i_, scaling roughly with the perturbative 
cross section. Figure 3 shows the logitudinal momentum spectrum integrated over all trans- 
verse momenta, and likewise Coulomb corrections persist to the highest and lowest values 
of q z . 

Given the decrease of Coulomb corrections with increasing beam energy one might ask, 
"At what 7 of colliding Pb beams would Coulomb corrections be relatively unimportant, 
say, less than 1% for the total cross section?" If for the purposes of reductio ad absurdum 
one takes the Lee Milstein formula as a reasonable order of magnitude approximation, then 
the answer is 7 = 10 43 . The point is that for any conceivable accelerator beyond LHC the 
Coulomb corrections to e + e~ pair production will still be significant. 
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FIG. 2: Positron transverse momentum spectrum for Au + Au at RHIC with 7 = 100. The filled 
circles are the exact calculation and the stars the perturbation theory. 

One can calculate momentum spectra to compare with the CERN SPS data. Since 
the CERN data comprise positrons uncorrelated with electrons, comparison with a full 
calculation of the positron momentum spectrum da(q) is appropriate. Figure 4 shows the 
data for a Pb projectile on a Au target. On the whole the perturbation theory curve (dashed 
line) perhaps seems closer to the data than to the solid full exact calculation. 

Figure 5 show an analagous comparison for a S projectile on a Au target. Again, the 
perturbation theory curve seems closer to the data, represented by the dot-dashed line. 
Figures 4 and 5 provide an illustration of the statement of the experimental authors, that 
the cross sections follow perturbative scaling. However, especially given the difficulty of the 
SPS experiment as described by the authors, the apparent lack of Coulomb corrections seen 
here needs to be verified in other ultrarelativistic heavy ion experiments. 
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FIG. 3: As in Figure 2 but for the positron longitudinal momentum spectrum for Au + Au at 
RHIC with 7 = 100. 

The first experimental observation of e + e~ pairs at RHIC has been published by 
STArJsi]]. Events were recorded where pairs were accompanied by nuclear dissociation. 
Comparison with perturbative QED calculations allowed a limit to be set "on higher-order 
corrections to the cross section, — 0.5ctqed < Act < §.2ijq ED at a 90% confidence level." 

The present technology of the properly regualarized exact computer code does not include 
an impact parameter representation and thus does not allow for evaluating a cross section 
where pair production is in coincidence with nuclear dissociation. Furthermore, the retarded 
propagators are not strictly appropriate when the range of both electrons and positrons 
are restricted such as the STAR data. However, a comparison of calculations in the STAR 
acceptance without nuclear dissociation is of interest as an indication of the relative difference 
between perturbation theory and the regularized exact result. In the STAR acceptance the 

13 




FIG. 4: Calculated positron momentum spectrum compared with the CERN SPS data for Pb + 
Au. The solid line is the exact calculation and the dashed line perturbation theory. 

exact result is calculated to be 17% lower than perturbation theory (as is coincidently true 
for the total RHIC e+e- cross section in Table I). This rough estimate Act = — 0.17ctqed is 
not excluded by the above STAR limit. 



IV. SUMMARY AND CONCLUSIONS 

A full numerical evaluation of the "exact" semi-classical total cross section for e + e~ 
production with gold or lead ions shows reductions from perturbation theory of 28% (SPS), 
17% (RHIC), and 11%(LHC). 

For large Z no final momentum region was found in which there was no reduction or an 
insignificant reduction of the exact cross section from the perturbative cross section. 
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FIG. 5: As in Figure 4 but with a S projectile. The dot-dashed line follows the authors' represen- 
tation of the CERN SPS data. 



The CERN SPS data cover a large part of the momentum range of produced positrons, 
and the present theory predicts a reduction of cross section at high Z from the perturbative 
result. That the CERN SPS data apparently do not show a reduction from perturbation 
theory is a puzzle. It would be of great interest to obtain more precise data on variation 
of heavy ion pair production cross sections with ion charge at RHIC or LHC. If the present 
apparent lack of evidence for Coulomb corrections in ultrarelativistic heavy ion e + e~ pair 
production were to be reproduced in other experiments it would provide a unique challenge 
to our theoretical understanding of strong field QED. 
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